Sympathetic cooling and growth of a Bose-Einstein condensate 
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We study two sets of rate equations for sympathetic cooling of harmonically trapped Bose gases. 
Calculations for mixtures of Na-Rb and Li-Cs show that both sets yield similar results for the cooling 
times. The equilibration rates are in fair agreement with each other and differ considerably from 
classical rates. The onset of Bose-Einstein condensation is rather sudden and non-exponential in 
time, and the growth of the condensate differs for the two different mixtures we studied. 
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I. INTRODUCTION 



The cooling of atoms in traps is an important tool 
in the study of ultracold Bose or Fermi gases. Usually, 
the last step in the cooling process is evaporative cool- 
ing. This technique, however, can only be used for suf- 
ficiently strongly interacting gases since thermodynamic 
equilibrium of the cooled gas must be (nearly) attained 
during the entire cooling process. There is a number of 
gases in which the interaction is too weak for evapora- 
tive cooling to work. Moreover, because of the exclusion 
principle the method fails at low temperature for identi- 
cal Fermions with parallel spins §. In such cases, one 
resorts to "sympathetic cooling": Another gas is cooled 
evaporatively and acts as the cooling agent for the gas to 
be cooled. 

The method of sympathetic cooling was proposed 
almost twenty years agofi || and has since found 
widespread application (see Refs. j|, ||). Exciting recent 
applications include the sympathetic cooling of atoms of 
different species: 41 K and 87 Rb §, 7 Li and 144 Cs §, of 
different isotopes of the same species: the bosonic species 
85 Rb cooled by 87 Rb ||] and 6 Li Fermions in a bath of 7 Li 
Bosons H [lCj , and the production of dual Bose-Einstein 
condensates (BEC) with sympathetic cooling|ll|, p^[ . 
Only a few theoretical descriptions of sympathetic cool- 
ing exist. The simplest descriptions are purely classi- 
cal while more sophisticated employ the quantum Boltz- 
mann equation or master equations. Delannoy et al. Jl^ ] 
and Mosk et al. Jli| derived a classical formula for the 
time dependence of the temperature and a classical for- 
mula for the cooling rate. Geist et al. have formu- 
lated sympathetic cooling of Fermi gases in terms of a 
quantum Boltzmann equation. The tremendous simpli- 
fication of the collision matrix elements achieved in this 
way precludes, however, a quantitative description of the 
cooling process. Lewenstein et al. derived a master 
equation for sympathetic cooling. Unfortunately, that 
master equation is too complex to be practically use- 
ful. Recently, the present authors simplified this master 
equation by means of decoherence and ergodicity argu- 
ments lla] and obtained two different sets of rate equa- 



tions that both describe sympathetic cooling of Bose or 
Fermi gases. These rate equations are simple enough to 
permit practical calculations while they contain detailed 
information about the quantum scattering processes be- 
tween the cooled gas and the bath. 

The purpose of this article is twofold. First, we want 
to compare the solutions of our rate equations with the 
much simpler expression for the classical cooling rate ob- 
tained by Delannoy et al. |Q and Mosk al. . Second, 
we use the rate equations to study the growth of a Bose- 
Einstein condensate for a system of 23 Na cooled by 87 Rb, 
and for 7 Li in a bath of 133 Cs. 

This article is organized as follows. In Section II we 
present the rate equations that we use. The results, i.e. 
the comparison between quantum and classical cooling 
rates and the studyof the growth of the condensate, are 
presented in Sect. III. Finally we give our summary. 



II. 



RATE EQUATIONS FOR SYMPATHETIC 
COOLING 



In this Section we present two different sets of rate 
equations that are derived in ref. Both sets of 

rate equations have been derived from a master equa- 
tion [la ] for sympathetic cooling by means of decoherence 
and ergodicity arguments. It is assumed that the cool- 
ing agent has a considerably larger mass than the gas 
to be cooled. The latter consists of Na (very weakly 
interacting) bosons in a spherically harmonic trap of 
frequency v. Single-particle trap orbitals have energies 
£m = hv(m x + m y + m z ) where rfi = (m Xl m yi m z ) de- 
notes the set of Cartesian quantum numbers, and the 
single-particle ground state energy is set to zero. 

The first set of rate equations describes how the mean 
occupation number of orbital too changes due to 
interactions with the cooling agent: 

±n Ao (t) = 2 £ r*^nrt(t)(n* (t) + l) 

m(^m ) 

- 2 rK(M*) + l)n*„(t).(l) 
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The form of these equations is intuitively obvious. Note 
that the rate equations (0) are of mean- field type and ne- 
glect correlations between occupancies of different single- 
particle orbitals. Particle number conservation implies 
that this approximation may fail when the occupancy of 
the ground state to = approaches the total number of 
particles in the gas to be cooled. Particular solutions 
to these mean-field type rate equations are discussed in 
Refs. @. 

The rate coefficients T™^ are input to the rate equa- 
tions (Q). For the computation of the rate coefficients 
r we assume that the bath particles are Boltzmann- 
distributed. A practical recipe for their computation is 
given in Ref.Q. We express the T's in units of 

u) = (327T 4 )- 1 A 3 B n B {a/l ) 2 [{M + to) 2 /{Mm)]v. (2) 

Here n B is the density of bath particles, A^ = 
(2nh 2 Pb/M) 1 / 2 their thermal de-Broglie wave length, 
Pb the inverse temperature, lo the oscillator length, and 
M and m are the mass of the cooling agent and the gas 
to be cooled, respectively. The interaction between bath 
and system is described in terms of the s-wave scattering 
length a. 

The second set of rate equations is based on a mi- 
crocanonical approach and describes how the probability 
Pm of having the gas to be cooled at total energy Mhv 
changes with time: 

dpM 



At. 



' / j n,rn 



x yPM+ a (nn[n,n + 1])a/ +q -pM<«ri["-m + 1])mJ (3) 

Here, afiv = £a~ £ m is the energy transfer, is the oc- 
cupancy of single-particle orbital to, and (n^ [n^ + 1]}m 
denotes the mean value of n^fn^ + 1] taken over the 
many-body states with fixed energy Mhv. They are thus 
microcanonical averages of occupancies. It is straight- 
forward to check that ^ M dpM(i)/di = 0. In Eq. (||), 
the pm's are the unknowns while the rate coefficients 
r and the quantities (n^n^ + 1])m are input quanti- 
ties defined in the framework of our model. The rate 
coefficients are identical to those used in the mean- 
field type rate equations (|l|). To compute the quanti- 
ties (ritflntft + 1]}m we neglect correlations by putting 
(nn[n,n + 1])m ~ («■«) A/K^m) m + 1] and replace the mi- 
crocanonical average by the grand canonical average 



n r7l {z,T) 

z{E) exp [— £js/ kT(E)] 



1 - z{E) exp [s^/nTiE)] 



(4) 



E=Mhv 



Expressions for the energy-dependent fugacity z(E) and 
temperature T(E) can be obtained numerically from 
the grand canonical partition function of non-interacting 
bosons in a harmonic trap. For details, we refer the 
reader to Ref. Jl6| . Note that the the microcanonical rate 
equations are expected to be valid also in the regime of 
BEC. 



III. RESULTS 

We study the cooling process for different tempera- 
tures of the bath and for different combinations of bath 
and system. The system always consists of Na = 400 
Bosons and the harmonic trap is cut off beyond the 
single-particle orbitals with energy Khv, K — 21. In 
the microcanonical approach, the homogeneous system 
of linear rate equations (||) can be put into matrix form. 
For stability reasons, the real parts of the eigenvalues of 
this matrix must be zero or negative. The equilibrium 
solution is determined by the eigenvalue zero. The equi- 
libration rate is given by the modulus of the real part 
of the eigenvalue with largest negative real part. The 
leading eigenvalues are computed using the sparse ma- 
trix solver Arpack |ll| . Within our numerical accuracy 
we found one zero eigenvalue and no eigenvalues with 
positive real parts. The results for the equilibration rate 
7 C q are listed in Table |[ The temperatures are chosen 
such that kT — Thv is close to the critical temperature 
while kT = 3hv yields a BEC at the end of the cooling 
process (k denotes Boltzmann's constant). 

We turn to the rate equations ([!]). We took the Bosons 
of the system A to be initially equally distributed over 
the degenerate single-particle orbitals with energy Khv, 
all other orbitals being empty. This initial configura- 
tion has the highest energy that can be realized in the 
finite trap we are considering. While details of the initial 
conditions are reflected in the short-time evolution they 
cannot affect the long-time behavior and the equilibra- 
tion rate. The rate equations (|l]) need to be modified in 
the regime of BEC; we therefore restrict the computation 
to the case kT = Ihv. We solve the rate equations nu- 
merically and extract the equilibration rate 7 eq from the 
solution at large times. The results are listed in Table Q. 
Note that the equilibration rates 7 oq and 7 oq depend only 
very weakly on the composition of system and bath for 
kT = Ihv. Comparing the rate 7 eq obtained in the micro- 
canonical approach with the rate 7 eq from factorization, 
we note that the former is about a factor w 1.7 larger 
than the latter. The time evolution however shows that 
the loss of energy at short times is practically identical 
in both approaches. Most of the energy is removed from 
system A during the early part of the cooling process. 
Therefore, both approaches yield comparable predictions 
for the cooling time. 

It is difficult to compare our full time dependence with 
that of the classical theory [jl2| 13 since we assumed a 
finite trap and, thus, different initial conditions. It is 
however meaningful to compare equilibration rates since 
these do not depend on initial conditions. In order to 
compare the equilibration rates, we have adjusted the 
formulas of Ref.|l2|, |l3| to the case of a homogeneous 
bath. In Table [jj, the classical rate 7 oq (which decreases 
quadratically with decreasing temperature) is seen to be 
about a factor 2 larger than the quantum-mechanical 
rate 7 eq ) below the condensation temperature, and a fac- 
tor 1.7 to 3 smaller in the regime of BEC. This casts some 
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Bath - System 


kT/Hu 


7c q /10 4 u; 


7eq/10 4 W 


7cq/10 4 W 


SY Rb - a Na 


7.0 


2.7 


1.6 


5.84 


8 'Rb - ^Na 


3.0 


3.1 




1.07 


i33 Cs - v Li 


7.0 


2.5 


1.7 


4.97 


i33 Cs - v Li 


3.0 


1.66 




0.913 



TABLE I: 

Comparison of quantum-mechanical equilibration rates 7 oq (obtained from the solution of Eq.([|)) and 7 cq (obtained from the 
solution of Eq. ([!])) with the classical rate 7 cq for several systems and temperatures. The critical temperature for onset of 
BEC is just below kT jTiv = 7. The classical equilibration rates differ considerably from the quantum mechanical rates. 



doubt on the possibility to describe sympathetic cooling 
in terms of the much simpler classical theory. 

It is particularly interesting to study the growth of a 
BEC. As pointed out above, the rate equations (Q) can- 
not be used in this regime. We chose lnT — ZTiv and 
followed the time evolution of the total energy E(t) = 
^"Em-^PmW by solving the rate equations (^) nu- 
merically. As initial conditions we took pkn a = 1 an d 
Pm = for M < KNa- This is the configuration with 
highest energy that might be realized in the open trap. 
Using the grand canonical partition function, we obtain 
the time-dependent temperature T(E) for the system, 
and the time-dependent ground state occupation number 
uq from Eq.(Q). Fig. [j] shows the result for two different 
systems. Note that the BEC grows faster in the Rb-Na 
system, which also has the larger equilibration rate [in 
units of u>, see eq. @]. The reason for this behavior is 
not easy to understand: Table [i] shows that equilibra- 
tion rates for both systems are quite similar at temper- 
ature kT ' /Tiv = 7 while they differ considerably below 
the critical temperature. Note that a naive (and classi- 
cal) explanation would link the difference in rates to the 
smaller momentum transfer during Li-Cs collisions and 
thereby to the larger mass difference. Our quantum me- 
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FIG. 1: Relative ground-state occupation as function of time 
for Li in Cs bath (full line) and Na in Rb bath (dashed line) 
at temperature kT = 3hv. 



chanical results show that the situation is more complex. 
The onset of BEC is rather sudden and non-exponential 
while the long time evolution is exponential. Qualita- 
tively similar observations have been made in numerical 
simulations of evaporative cooling (l^, |2^, |2l[] . One does 
not obtain 100% BEC since the final temperature is only 
about half the condensation temperature and the system 
is not of macroscopic size. 



IV. SUMMARY 

In summary, we have used two sets of rate equations 
to describe sympathetic cooling of Bose gases by a much 
heavier gas. Practical calculations in the systems 23 Na- 
87 Rb and 7 Li- 133 Cs show that the cooling times from 
both rate equations agree well with each other. Our 
quantum-mechanical equilibration rates are in fair agree- 
ment with each other but differ considerably from the 
classical rates. This is unfortunate since the classical ap- 
proach is much simpler than the quantum-mechanical de- 
scription presented in this work. We studied the growth 
of the Bose-Einstein condensate and found that the con- 
densate in the system 23 Na- 87 Rb grows faster than in 
7 Li- 133 Cs. At the onset of Bosc-Einstcin condensation 
we observe a rather sudden and non-exponential growth 
of the condensate. 



Acknowledgments 

The authors are grateful to A. Mosk and M. Wei- 
demiiller for suggesting this problem, and for many stim- 
ulating discussions. They also thank S. J. Wang for valu- 
able suggestions. T.P. thanks the Max-Planck-Institut 
fur Kernphysik, Heidelberg, for its hospitality during 
the initial stages of this work. Oak Ridge National 
Laboratory is managed by UT-Battelle, LLC for the 
U.S. Department of Energy under Contract DE-AC05- 
00OR22725. 



4 



[15] 



B. DeMarco et al, Phys. Rev. Lett. 82 (1999) 4208. 

D. J. Wineland, R. E. Drullinger, and F. L. Walls, Phys. 
Rev. Lett. 40 (1978) 1639. 

R. E. Drullinger, D. J. Wineland, and J. C. Berquist, 
Appl. Phys. 22 (1980) 365. 

D. J. Larson et al, Phys. Rev. Lett. 57 (1986) 70. 

Atomic Physics 9, edited by R. S. van Dyck, Jr., and E. 

N. Fortson, World Scientific, Singapore 1985. 

G. Modugno et al, Scienc e 294, 1320 (200 1). 

M. Mudrich et al, arXiv: physics/0111213 . 

I. Bloch et al., Phys. Rev. A 64, 021402 (2001). 

A. G. Truscott et al, Science 291 (2001) 2570. 

F. Schreck et al., Phys. Rev. A 64 (2001) 011402(R). 

C. J. Myatt et al., Phys. Rev. Lett. 78 (1997) 586. 

G. Delannoy et al, Phys. Rev. A 63 (2001) 051602(R). 
A. Mosk et al., Appl. Phys. B 73, 791 (2001). 

W. Geist, L. You, and T. A. B. Kennedy, Phys. Rev. A 
59 (1999) 1500. 

M. Lewenstein, J. I. Cirac, and P. Zoller, Phys. Rev. A 



51 (1995) 4617. 

[16] T. Papenbrock, A. N. Salgueiro, and H. A. Weidenrmiller, 
Phys. Rev. A 65, 043601 (2002). 

[17] S. J. Wan g M. C. Nemes, A N. Salgueiro, H. A . Wei- 
denmiiller, |cond-mat/0110293| ; |cond-mat/01 11413 . 

[18] R. B. Lehoucq, D. C. Sorensen, and C. Yang, ARPACK 
User's Guide (SIAM 1998) ISBN 0-89871-407-9. 

[19] O. J. Luiten, M. W. Reynolds, J. T. M. Walraven, Phys. 
Rev. A 53, 381 (1996). 

[20] M. Holland, J. Williams, and J. Cooper, Phys. Rev. A 
55, 3670 (1997). 

[21] C. W. Gardiner et al, Phys. Rev. Lett. 79, 1793 (1997); 
C. W. Gardiner and P. Zoller, Phys. Rev. A 58, 536 
(1998); C. W. Gardiner et al, Phys. Rev. Lett. 81, 5266 
(1998); M. D. Lee and C. W. Gardiner, Phys. Rev. A 
62, 033606 (2000); M. Kohl et al, Phys. Rev. Lett. 88, 
080402 (2002). 



